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Abstract 

A selective chronological survey of Darboux transformations as related to 
supersymmetric quantum mechanics, intertwining operators and inverse 
scattering techniques is presented. Short comments are appended to each 
quotation and basic concepts are explained in order to provide a useful 
primer. 



Contents 

Chronological set 1: DTs (as covariance of Sturm-Liouville problems) 

Chronological set 2: Intertwining (transformation) operators (the more general 
operator language) 

Chronological set 3: Related inverse quantum scattering approaches 
Some abbreviations: 

1. Darboux transformation(s) - DT(s) 2. Factorization method - FM 

3. Inverse Quantum Scattering - IQS 4. One-dimensional - ID 

5. Quantum mechanics - QM 6. Schroedinger equation - SE 

7. Supersymmetry, supersymmetric - SUSY 



Possible combinations are allowed. 
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CHRONOLOGICAL SET 1 
Darboux Transformations 
(covariance of Sturm-Liouville problems) 

1882: Sur une proposition relative aux equations lineaires 

This proposition, which history proved to be a notable theorem, is provid ed by 
G. Darboux in Compt. Rend. Acad. Sci. 94, 1456 ( |physics/9908003| ), and 
is also to be found in his course "Theorie des Surfaces", II, p. 210 (Gauthier- 
Villars, 1889). Darboux's result was not used and recognized as important for 
a long time, and in fact has been only mentioned as an exercise in 1926 by Ince. 



1926: Ordinary Differential Equations (Dover) 



In his famous Dover book, E.L. Ince publishes Darboux's theorem in a section 
of "Miscellaneous Examples" at page 132, together with other two particular 
examples of the theorem related to the free particle and Poschl- Teller potential 
(the latter belonging to Ince), Exercises 5, 6 and 7, respectively. The exercises 
in Ince follow closely the formulation of Darboux (see §1). 

§1: Exercise 5 at page 132 of Ince's textbook 
Prove that if the general solution u = u(x) of the equation 



— = (0(x) + h)u 



is known for all values of h, and that any particular solution for the particular 
value h = hi is u = f(x), then the general solution of the equation 



for h =/= hi is 



y = u (x) - u(x)l—- 



1940-1941: The factorization of the hypergeometric equation 



In Proc. Roy. Irish Acad. 47 A, 53 (1941) flphysics/9910003D , E. Schrodinger 
factorizes the hypergeometric equation, finding that there are several ways of 
factorizing it. This was a byproduct of his FM (published during 1940-1941 in 
the same review) originating "from a, virtually, well-known treatment of the 
oscillator", i.e., an approach that can be traced back to Dirac's creation and 
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annihilation operators for the harmonic oscillator and to older factorization 
ideas in the mathematical literature. 

1951: The FM 

L. Infeld and T.E. Hull present the classification of their factorizations of linear 
second order differential equations in Rev. Mod. Phys. 23, 21. 

1955: Associated Sturm- Liouville systems 

M.M. Crum publishes an important iterative generalization of Darboux's result 
in Quart. J. Math. 6, 121 ( |physics/9908019| ), without any mention of Darboux 
(see §2). 

§2: Cram's iteration 

Let ipi,ipa, •••ipff be solutions of a given Schrodinger equation —D 2 ipi + uipi = 
Xijji, for fixed, arbitrary constants A = A 1; X 2 , ...Ajy, respectively. Define the 
Wronskian determinant W of k functions /i, f%, fa by 

W(fi,...,f 2 ) = detA, A l3 = ^-J±, i,j 1.2 k. 

Then, the transformations (Cram's formulas) 

u[N] = u-2D 2 lnW(ip ll ... 1 ip N ) 
are covariant ones, i.e., one can write the following SL equation: 

-D 2 ^[N] + u[N]ip[N] = \ip[N] . 

Darboux's result of 1882 may be seen as the case N = 1. In other words, Crum 
presented the succesive (iterative) DTs in compact formulas. 

1957: Krein's approach to generalized DT 

M.G. Krein publishes extensions of Crum's iteration in D.A.N. SSSR 113, 970, 
using a different perspective and with more elaborations. In particular he pro- 
vides an important theorem on the sign of the Wronskian. 



1979: Concept of DT 
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V.B. Matveev uses for the first time the expression DT (directly in the title) in 
his series of papers in Lett. Math. Phys. 3, 217, 425, in which he extended the 
range of applicability of the concept to differential-difference and difference- 
difference evolution equations. 

1981: Dynamical breaking of SUSY 

In section 6, "Some models", of his renowned paper Nucl. Phys. B 185, 513, 
E. Witten introduces SUSY QM as a toy model for supersymmetry breaking in 
quantum field theories. As a matter of fact, the whole SUSY breaking paradigm 
is nothing but a classification of vacuum states based on the SUSY QM toy 
model. The SUSY breaking is considered as a sort of "phase transition" with 
the order parameter being the Witten index, defined as the grading operator 
t = (— 1) W/ , where Nf is the fermion number operator. For the ID SUSY 
QM Witten's index is the third Pauli matrix 03, which is +1 for the bosonic 
sector and -1 for the fermionic sector of the ID quantum problem at hand, ft 
is also quite common to call a particular Riccati solution as a (Witten) "su- 
perpotential" . Papers that now are standard references are published during 
1982-1984. It is known that H. Nicolai was the first to write the simple SUSY 
QM matrix algebra in 1976, whereas a relativistic graded algebra has been first 
written by Yu.A. Gol'fand and E.P. Likhtman in JETP Lett. 13, 323 (1971) 
and since then supersymmetric gauge theories have been main stream research. 

1981: DTs and nonlinear evolution equations 

V.B. Matveev and M. Salle apply for the first time DT to nonlinear equations 
in D.A.N. SSSR 261, 533. At present, there are hundreds of papers in the area. 

§3: DTs and Backlund Transformations 

A well-know transformation in nonlinear physics is the Backlund transforma- 
tion (BT) [A.V. Backlund, Math. Ann. 9 (1876), 297; 19 (1882), 387], that 
has been first applied to sine-Gordon equation. When BT operators are ap- 
plied to N soliton solutions one gets N+l soliton solutions. In the Korteweg-de 
Vries case, DTs may be seen as BTs for instantaneous solitons, since instanta- 
neous KdV sech 2 solitons may be considered at the same time as Schrodinger 
potentials of a well-known SUSY QM problem. 

On the other hand, in the Lax representation of KdV equation, the first Lax 
operator is of Schrocdinger type, a fact showing the importance of Darboux 
covariance in nonlinear physics. 

1983: Derivation of exact spectra of the SE by means of SUSY 
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In JETP Lett. 38, 356, L.E. Gcndenshtein introduces the important concept 
of shape invariance (SI) in SUSY QM. 

§4: The SI property 

SI is a property of some classes of potentials with respect to their parameter(s), 
say a, and reads 

V n+1 (x,a n ) = V n (x,a n+1 ) + R(a n ) 

where R is a remainder. This property assures a fully algebraic scheme for the 
spectrum and wavefunctions. Fixing Eq = 0, the excited spectrum is given by 

n+l 

and the wavefunctions are obtained from 

n 

V>„(a;,cn) = Y[ A + (x,a k )ip (x,a n+1 ) 
k=i 



1984: FM and new potentials with the oscillator spectrum 

In J. Math. Phys. 25, 3387, B. Miclnik provides the first application of the 
general Riccati solution to the harmonic oscillator, noticing also the similar- 
ity to the Abraham-Moses class of isospectral potentials in the area of inverse 
scattering. In the same year, D. Fernandez gives a second application to the 
atomic hydrogen spectrum in Lett. Math. Phys. 8, 337, whereas M.M. Nieto 
clarifies further the inverse scattering aspect of Miclnik's construction in Phys. 
Lett. B 145, 208. The procedure may be seen as a double Darboux transfor- 
mation in which the general Riccati (superpotential) solution is involved (see 
§7). Thereby it will be denoted by the acronym DDGR (double Darboux gen- 
eral Riccati). 

1984-1985: FM and DTs for multidimensional Hamiltonians 

In papers published in Teor. Mat. Fiz. 61 and subsequent Phys. Lett. A 
versions, A. A. Andrianov, N.V. Borisov and M.V. Ioffe discover the relation 
between SUSY QM and DTs while playing with matrix Hamiltonians in SUSY 
QM. 

1985: Exactness of semiclassical bound state energies for SUSY QM 
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In Phys. Lett. B 150, 159, A. Comtet, A. Bandrauk and D.K. Campbell were 
the first to study the WKB features of SUSY QM introducing a SUSY QM 
WKB formula QE -W 2 (y)] 1 / 2 dy — nwh, where W is Witten's superpotential 
and a and b are turning points. This research line has a rich publishing output. 

1988: PARASUSY QM 

V.A. Rubakov and V.P. Spiridonov introduce PARASUSY QM in Mod. Phys. 
Lett. A 3, 1337, involving supercharges of order-three nilpotency. Later, J. 
Beckers and N. Debergh found a different algebra of the same type [Nucl. 
Phys. B 340, 767 (1990)]. There are many generalizations, making this ten- 
year topic quite active and interesting from the application point of view. 

1991: DTs and Solitons (Springer) 

V.B. Matveev and M. Salle publish the first (excellent) book (112 pp, 197 refs.) 
focusing on DTs and their relation with soliton (mathematical) physics. The 
DTs are defined as covariant properties of the SE. Many types of DTs (not 
mentioned here) are presented in this book in a concise manner. 

§5: Darboux Covariance 

The statement of the Darboux theorem can be interpreted as the Darboux 
covariance of a Sturm-Liouville equation 

-Ipxx +U1p = \lp 

by which one should understand that the following DT 

tp —> tp[i] = (D - <ji)v> = Vx - ctiV 1 = 
i>i 

u — > u[l] = u — 2u\ x =u — 2D 2 ln^i 

where o\ — ipixip" 1 , i-e. it is the logarithmic derivative, and W is the Wron- 
skian determinant passes the SL equation to the (Darboux isospectral) form 

-V*x[l]+«[1M1] = 

When DTs are applied iteratively one gets Crum's result. One can also say 
that the two SL equations are related by a DT. 
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1994: A modification of Crum's method 

In Theor. Math. Phys. 101, 1381, V.E. Adler discusses a double Darboux 
transform for consecutive eigenfunctions i\>k and V'fe+i as transformation func- 
tions, producing a nonsingular potential missing the two levels E = k and 
E = k + 1 in its discrete spectrum (thus anharmonic). According to Samsonov 
this result is implicit in Krein's theorem on the sign of the Wronskian. Very 
recently, D.J. Fernandez, V. Hussin and B. Mielnik [Phys. Lett. A 244, 309, 
(1998)] gave an interesting combination of Adler's method and DDGR. 

1994: Coherent states for isospectral oscillator Hamiltonians 

In J. Phys. A 27, 3547, D.J. Fernandez, V. Hussin and L.M. Nieto provide the 
first discussion of coherent states for Darboux transformed systems. In this im- 
portant topic there are recent significant results due to Bagrov and Samsonov 
[J. Phys. A 29, 1011 (1996)] and Samsonov [J. Math. Phys. 39, 967 (1998)]. 
In the latter paper, Samsonov showed that in some cases the distorsion of the 
phase space due to DT is calculable. 

1995: SUSY and QM 

In Phys. Rep. 251, 267, F. Cooper, A. Khare and U. Sukhatme publish the 
latest (at this time) review on SUSY QM. Even though this is a comprehen- 
sive work with 265 references, they cite 16 omitted topics and are forced to 
accept the following: "So much work has been done in the area of SUSY QM 
in the last 12 years that it is almost impossible to cover all the topics in such 
a review" . For comparison, what may be considered (in a limited sense) as the 
first SUSY QM review paper written by L.E. Gendenshtein and I.V. Krive in 
1985 has 65 references in various areas. 

1996: SUSY Methods in Quantum and Statistical Physics (Springer) 

G. Junker publishes his Habilitazion at Erlangen as a book (124 pp) with more 
than 300 references. 

1996: DTs for time- dependent SE (TDSE) 

In Phys. Lett. A 210, 60, V.G. Bagrov and B.F. Samsonov provide an im- 
portant extension of DTs to nonstationary SEs by means of intertwining. A 
direct (less general) factorization approach for the time-dependent Pauli equa- 
tion has been given by V.A. Kostelecky, V.I. Man'ko, MM. Nieto, and D. 
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Truax, in Phys. Rev. A 48, 951 (1993); see also, V.M. Tkachuk, J. Phys. A 
31, 1859 (1998). 

1997: DT for Dime eqs with (1+1) potentials 

A.V. Yurov studies DT for Dirac equations in Phys. Lett. A 225, 51. 
1997: DT of the SE 

In Phys. Part. Nucl. 28, 374 (1997) [Fiz. Elem. Chastits At. Yadra 28, 951 
(1997)] V.G. Bagrov and B.F. Samsonov publish a first review of the results in 
QM obtained by means of DT. There are 83 references. 

§6: DTs and SUSY QM (Matveev and Salle) 

Proposition: Witten's SUSY QM is equivalent to a single DT. 

Proof: Consider two Schrodinger equations 

-L> 2 ?A + uip = Xip 

-D 2 <f> + v<t> = \<f> 
related by DT, i.e., v — u[l] and <j> = tp[l]. 

Notice now that the function <j>i = i})^ 1 satisfies the Darboux-transformed 
equation for A = Ai. 

If now one uses the second (transformed) equation as initial one and perform 
the DT with the generating function <j>i, one just go back to the initial u 
equation. Thus, one can define an inverse DT, that follows in a clear way from 
the direct one: 

u = v - 2D 2 \n<p 1 = v[-l] = v- 2D 2 ln?/>r 1 

i> = (<fe - (Ai - A) = (V + t±<^ (Ai - A) 

If the sigma notation for the logarithmic derivative is introduced, i.e., 

_ ?hx _ <hx 

the Riccati (SUSY QM) representation of the Darboux pair of Schrodinger 
potentials is obtained 

u = v[-l] = a x + a 1 + Ai 
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v = u[l] = —o x + a + Ai 

It is now easy to enter the SUSY QM concept of supercharge operators. For 
that, one employs the factorization operators 

B+ = -D + a, B~ = D + a 

They effect the wavefunction part of the direct and inverse DT, respectively. 
Moreover 

B+B- = -D 2 + v-\ 1 

B-B+ = -D 2 + u-\ 1 

Thus, the commutator [B + ,B~] = v — u = — 2D 2 ln?/'i gives the Darboux 
difference in the shape of the Darboux-related potentials. Introducing the 
Hamiltonian operators 

H+ = B-B+ + X 1 
H~ = B+B~ +Ai 

one can also interpret the B operators as factorization ones and write the 
famous matrix representation of SUSY QM, as well as the simplest possible 
superalgebra. 

The factorizing operators in matrix representation are called supercharges in 
SUSY QM, and are nilpotent operators 



Q = A-(t + 





A- 



and 



« + - ^- - ( i - 4 ; ) 



er_ = y q J J and a + = y ^ jj j arc Pauli matrices. In this realization, 
the matrix form of the Hamiltonian operator reads 



H 



A+A- \ _ / H- 
A-A+ ) ~ \ H+ 



defining the partner Hamiltonians as diagonal elements of the matrix one. They 
are partners in the sense that they are isospectral, apart from the ground state 
4> gr - of H_, which is not included in the spectrum of H + . 

Remark: Denoting A\(x) = — 2D 2 In ipi{x) + Ai, one can see that for a Gaus- 
sian of width a 2 one gets Ai = 2cr _2 + Ai, and in the case of Gaussians of 
constant width (harmonic ground states) Ai is just a constant shift of the fac- 
torization constant. 
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§7: The DDGR Construction (Rosu) 



DDGR offers an interesting possibility to construct families of potentials strictly 
isospectral with respect to the initial (bosonic) one, if one asks for the most 
general superpotential (i.e., the general Riccati solution) such that V+(x) = 
w g + "a?' wncrc V + is the fermionic partner potential. It is easy to see 
that one particular solution to this equation is w p = w(x), where w(x) is the 
common Witten superpotential. One is led to consider the following Riccati 
equation w 2 + = w 2 + , whose general solution can be written down as 
w g (x) = w p (x) + where v(x) is an unknown function. Using this ansatz, 
one obtains for the function v(x) the following Bernoulli equation 



dv(x) 
dx 



- 2v(x) w p (x) = 1, 



that has the solution 



v(x) 



X (x) + A 
ug(x) ' 



(1) 



(2) 



where Io( x ) = f* u o(y) dy, (c = — oo for full line problems and c = for half 
line problems, respectively), and A is an integration constant thereby considered 
as a free DDGR parameter. Thus, w g (x) can be written as follows 



w g (x; A) = w p (x) 



dx 



ln(2b(a:)+A) 



_d 

"dx 



w p (x) + cr (A) 
u (x) 



In 



,J (x) + A / 

Finally, one easily gets the V- (x; A) family of potentials 

dw g (x; A) 



V_(x;A) = w 2 (x;A)- 



dx 



(3a) 
(36) 
(3c) 

(4a) 



= V_(x)-2 



dx 2 
V_(x)- 



ln(J (x) + A) 
2ao,x(A) 



= V_(x)- 



4u (x)u (x) 2 U 4(x) 



+ 



J (x) + A (2o(x) + A) 2 



(46) 
(4c) 

(Ad) 



All V_ (x; A) have the same supersymmetric partner potential V+ (x) obtained 
by deleting the ground state. They are asymmetric double-well potentials that 
may be considered as a sort of intermediates between the bosonic potential 
V_(x) and the fermionic partner V + (x) = V_(x) — 2cr ,x(x). From Eq. (3c) 
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one can infer the ground state wave functions for the potentials V_(x;A) as 
follows 

u (x;A) = f(A) 7 (X ] (5) 

where f(A) is a normalization factor that can be shown to be of the form 
f(A) = y/ A(A + 1). One can now understand the double Darboux feature of the 
DDGR by writing the parametric family in terms of their unique "fermionic" 
partner 

V_ (x; A) = V+(x) - 2-^ In ( , (6) 



dx 2 \ uo (x; A) t 

which shows that the DDGR transformation is of the inverse Darboux type, 
allowing at the same time a two-step (double Darboux) interpretation, namely, 
in the first step one goes to the fermionic system and in the second step one 
returns to a deformed bosonic system. 



CHRONOLOGICAL SET 2 
Intertwining (Transformation) Operators 
(the more general operator language) 



1938: Intertwining discovered 



J. Delsarte introduces the notion of transformation (transmutation) operator 
in Sur certaines transformationes fonctionelles relative aux equations lineaires 
aux derivees partielles du 2nd ordre, Comp. Rend. Acad. Sci. (Paris) 206, 
1780 (physics/9909061), [see also, Sur une extension de la formule de Taylor, 
J. Math. Pures et Appl. 17, 213 (1938)], but it is only 18 years later that 
he with J.L. Lions expound more on the method. J. Delsarte, Colloq. Int. 
Nancy, p. 29 (1956); J.L. Lions, ibidem, p. 125 (1956); J. Delsarte, J.L. Lions, 
Comment. Math. Helv. 32, 113 (1956). 



§8 What is intertwining ? 



Two operators Lq and L± are said to be intertwined by an operator T if 

L\T = TL Q . (1) 

If the eigenfunctions ipo of Lq are known, then from the intertwining rela- 
tion one can show that the (unnormalized) eigenfunctions of L\ are given by 
ipi = TifQ. The main problem in the intertwining transformations is to con- 
struct the transformation operator T. ID QM is one of the simplest examples 
of intertwining relations since Witten's transformation operator T qm = T± is 
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just a first spatial derivative plus a diffcrcntiable coordinate function (the su- 
perpotential) that should be a logarithmic derivative of the true bosonic zero 
mode (if it exists), but of course higher-order transformation operators can be 
constructed without much difficulty. 

Thus, within the realm of ID QM, writing T\ = D — where u is a true 

bosonic zero mode, one can infer that the adjoint operator T\ = —D — ^- 
intertwines in the opposite direction, taking solutions of L\ to those of Lq 



In particular, for standard ID QM, L = H- and L\ = H + and although the 
true zero mode of H- is annihilated by Ti, the corresponding (unnormalizcd) 
eigenfunction of H+ can nevertheless be obtained by applying T\ to the other 
independent zero energy solution of H- . 

1973: Theory of generalized shift operators (Nauka) 

This is one of the remarkable books of B.M. Levitan edited by the publishing 
house Nauka. 

1978: Applications of a commutation formula 

P.A. Deift presents applications of the so-called "commutation formula" (that 
can be found, e.g., in the book of S. Sakai, C* -Algebras and W* -Algebras, 
Springer, 1971). In the paper of Deift in Duke Math. J. 45, 267 (1978), section 
4 contains the application of the commutation formula to ordinary differential 
operators. Many of the results in that section have been reproduced later in 
SUSY QM style, e.g. by Sukumar. 

1986-1987: Isometric operators, isospectral Hamiltonians, and SUSY QM 

In Phys. Rev. D 33, 2267, D 36, 1103, D.L. Pursey uses intertwining in his 
combined procedures of generating families of strictly isospectral Hamiltonians, 
starting from the Marchenko inverse scattering equation. 

1991: Intertwining of exactly solvable Dirac eqs. with ID potentials 

In Phys. Rev. A 43, 4602, A. Anderson applies matrix intertwining relations 
to the Dirac equation showing that their structure is described by an N— 4 su- 
peralgebra. 




1V1 • 



(2) 



1995-1998: Intertwining widely used 
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Intertwining is already well known to many active authors in SUSY QM, who 
are playing with higher-order generalizations. But, as always, the most impor- 
tant (at least for standard quantum mechanics) are the simplest cases, namely 
the Darboux first-order intertwining operators. 



§9: Using intertwining: DTs for TDSE (Bagrov and Samsonov) 



The DT for the TDSE are based on the intertwining relation 

T(id t - H ) = (id t - HJT , 

where 



H, 



-d£ + Vi(x,t), i = 0,l, 



and T is a first-order diff. transformation operator of the form 

T = L 1 (x,t)d x +L {x,t). 



(1) 
(2) 

(3) 



It follows immediately from the intertwining relation that if i/jq solves the TDSE 
with Hamiltonian Ho, then ipi = Tipo will solve the TDSE with Hamiltonian 
Hi. It is also easily verified that the intertwining relation will be satisfied if 
and only if 



T = Li{d x + xx) , Vi = V + 2 Xxx + iQogLOt 



(4) 



where e~ x is a solution of the TDSE with potential Vo, and Li = Li(t) is an 
arbitrary function. The transformed potential Vi(x, t) is a real- valued function 
if and only if 

ImXxxx = 



and 



\Li\ = cxp 



2 / lmx xx {x, s)ds 

J to 



Without loss of generality, one can assume that Li is real and positive, and is 
therefore given by the right-hand side of the above equation. 

Just as in the time-independent case, the DT for the TDSE can be inverted. 
Thus, if ipi is a solution of the TDSE with potential Vi given by Eq.(4), the 
function 



ipo{x,t) = 



>-x(x,t) 



i x ^My,t)dy + co{t) 



(5) 



with Co(t) given by 

[•t e x{xa,s) 

c (t)=iLi(t) / (il>i,x(xo,s) -Xx{x ,s)ipi(x ,s)) 

J to L i( s ) 



ds 
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is a solution of the T DSE with potential Vq. As remarked by F. Finkel et al. 
in math-ph/9809013 if the factor L\ is taken as unity, the mapping ip\ i— > 
given by Eq.(5) reduces to the non-local transformation considered by Bluman 
and Shtelen in J. Phys. A 29, 4473 (1996). 



CHRONOLOGICAL SET 3 
Related Inverse Quantum Scattering (IQS) Approaches 

In the case of classical dynamics, the inverse problem just means to determine 
the force acting on a macroscopic body from the features of the trajectory. In 
this sense, one of the oldest inverse problem has been Newton's problem of 
determining the force on the planets from the Kepler properties (laws) of their 
movement. 

In the realm of differential operators, an inverse spectral method generally 
means to determine the (linear) operator from some given spectral data. Roughly 
speaking, in the case of a Sturm-Liouville operator, this would reduce to get- 
ting the potential. Thus, there is a certain parallel to the classical case. As 
well known, in the IQS approach one uses integral transformation operators. 
The key point for the construction of the integral operators are the famous 
Gel'fand-Levitan equation and Marchenko equation. The IQS research started 
in the second half of the 1940s and the 1950s have been a real boom. 



1954-1956: Krein's approach to IQS 



M.G. Krein reports a new IQS approach in three D.A.N, notes, DAN SSSR 97, 
21 (1954); 105, 433 (1955); 111, 1167 (1956). 



§10: Krein's IQS 



As explained to us by Chadan and Sabatier in their book (1977, 1989), Krein's 
approach is the following way to solve the I = inverse scattering problem. 
The SE 

-D 2 y + V(r)y = k 2 y , 
where D = 4-, is substituted by the equivalent system 

Dy + A(r)y = kz 

—Dy + A(r)z = ky . 
The two functions V and A are connected by the Riccati equation 

V = -DA + A 2 . 
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From the analytical properties of the Schrodingcr regular solution 4>(k, r) and 
the Wiener-Paley theorem, one can obtain the following representation 



where A{r) = 2r 2r (2r). Moreover, for any fixed value of r, the function T 2r (t) 
is a solution of the Fredholm integral equation 



where, H{t) is related to the Jost function F(k) in the following way 



The scheme to solve the inverse scattering problem according to Krein is first 
to build the function H(t), then to get T2 r (t), leading immediately to A(r), 
and finally obtaining the potential V(r) from the Riccati equation. 

1980: Changes in potentials due to changes in the point spectrum: Anharmonic 
oscillators with exact solutions 

Using the Gel'fand-Levitan equation, P.B. Abraham and H.E. Moses introduce 
a procedure of deleting and adding bound states, which is "almost" equivalent 
to DDGR [Phys. Rev. A 22, 1333]. 

§11: Comparison between strictly isospectral techniques within IQS 

According to A. Khare and U. Sukhatme [Phys. Rev. A 40, 6185 (1989)], the 
difference between DDGR, Abraham-Moses, and Pursey schemes when used 
in the IQS context to generate one-parameter strictly isospectral potentials on 
the base of an initial zero mode tpo lies merely in the employed function used 
for the (inverse) steps of readding the zero- mode as follows (I = /_ iPo(y)dy): 

For DDGR, the function is <p = tp^ , being a node- free solution of the first step 
transformed equation. 

For Pursey 's method, it is v = tpo/ 1, such that v is a solution of the first step 
zero-energy SE that vanishes at +oo. 

For the Moses- Abraham procedure, one should use u = ipo/(l — I), where u is 
a solution of the first step zero-energy SE vanishing at — oo. 

One can obtain five independent families of strictly isospectral potentials when 
combining the three procedures. However, only the DDGR method leads to 
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reflection and transmission amplitudes identical to those of the original poten- 
tial, showing the complete degeneracy produced by such a construction. 

1985: First issue of Inverse Problems 

Birth year of the IOP review Inverse Problems. 

There are many good books in the IQS field, e.g., those of K. Chadan and 
P.C. Sabatier (1977, 1989), B.M. Levitan (1984), V.A. Marchcnko (1986), B.N. 
Zakhariev and A. A. Suzko (1985). According to Levitan, the generalized use 
of intertwining operators in IQS is due to the works of V.A. Marchenko. 

1985: SUSY QM and the Inverse Scattering Method 

In J. Phys. A 18, 2937, C.V. Sukumar provides an important contribution to 
the understanding of the connections between the two approaches. 

1993: Bound states in the continuum (BSICs) from SUSY QM 

In Phys. Rev. A 48, 3525, J. Pappademos, U. Sukhatme, A. Pagnamenta are 
the first to apply DDGR to get BSICs. They also study two-parameter families 
of SUSY BSIC potentials. 

1995: SUSY transformations of real potentials on the line 

In J. Phys. A 28, 5079, J.-M. Sparenberg and D. Baye present an exhaus- 
tive work on SUSY iterations showing the power of the method in generating 
isospectral potentials. 

1994: Exactly solvable models for the SE from generalized DTs 

In J. Phys. A 27, 2605, W.A. Schnizer and H. Lccb show in the IQS context 
that integral transformations with a degenerate kernel are equivalent to differ- 
ential ones. 

1995: On the equivalence of the integral and the differential exact solution gen- 
eration methods for the ID SE 

In J. Phys. A 28, 6989, B. Samsonov shows the same equivalence in a very 
transparent way, giving an integral form of Crum's iteration. 
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Conclusion 

A scientific review covering 116 years, even when the chronological style is used, 
cannot have any claim of completeness and people reading it may see a lot of 
missing material in general due to the natural biases acting on the author(s). 
In my case, a strong bias has been a dead line forcing mc to leave out many 
active authors and their results in a number of important topics not quoted 
above, such as the connections between partially algebraic quantum problems 
(quasi-exactly solvable ones) and SUSY QM, or intricate connections going 
directly to the core of more sophisticated theories. Nevertheless, even from 
this brief chronological collection one can get a partial grasp of what has been 
done and a feeling of what could be done in the future. 
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